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for every plate than to consider whether it will be necessary to 
do so in particular cases. 

* In the calculation of the tables from which the diagrams 
were drawn I have received valuable assistance from Miss Julia 
Bell and Miss Hilda Gibb, of Girton College, to whom my 
acknowledgments are due. 

We are indebted to the Government Grant Fund of the 
Royal Society for a grant to secure their assistance in the reduc¬ 
tion of the Eros photographs. 


Cambridge Observatory : 1903 January 3. 


Note on the Use of Mr. Aldis ’ Tables of the Function i (0 + cos 0 ) 
in Determining the Elements of an Orbit. By H. C. 
Plummer, M.A. 


1. In the second part of his explanation* of the use of 
tables of ±(d + cos 6 ) Mr. Aldis has shown that they can be 
applied to the problem of finding the elements of an elliptic orbit 
when two focal distances, the angle between them, and the time 
occupied in passing from one position to the other are known. 
An inspection of the example which he has worked out makes it 
evident that the practical utility of the method may be greatly 
enhanced if a convenient way can be devised of finding preliminary 
approximations to the values of the quantities required. It is 
the object of this note to supply the want. A graphical method 
will be first described ; but this is to be regarded chiefly as the 
basis of the analytical method of approximation which follows. 

2. By the use of Lambert’s theorem the problem is reduced 
to the solution for a, e, and £ of this system of three equations : 


— = (e — sin e) — (8— sin 8) 
a 7 

r 4* r' — c 


1 — cos 8 = 


I — COS £ = 


2a 
r 4- r' + c 
2a 


(1) 

... (2) 

- ( 3 ) 


The notation is that of Mr. Aldis’ paper. When a is known 
the eccentricity e can be easily deduced. It will be convenient 
to write 

o-j 2 = r+r' + c ; a 2 2 = r+r'—c. 

Now if 

x x = c— sin 2 ; = 1 — cos Z 

x 2 = e — sin e ; y z = 1 — cos e 

* Monthly Notices, vol. lxii. p. 638. 
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148 Mr. Plummer , Note on the Use of exiii. 3, 

x„ and x 2 , y 2 are the coordinates of points lying on the cycloid 
whose equations are # 

x = 6 — sin 0 ; y — 1 — cos 6 

The coordinates are also connected by the relations 

ht = o?(x 2 — xf) 

2 ay x = o - 2 2 ; 2 ay 2 = <r x 2 

If a be eliminated we have 

2 \/ 2ktyft = rrp>{x 2 — xf) 

2\f 2ktyf = cti 3 (^ 2 — x f) 

3. These equations have a geometrical interpretation. Let 
p and q (fig. 1) be the points (x t , yft) and (x 2 , y 2 ) on the cycloid 



ACB. The meaning of the equations is that semi-cubical 
parabolas with known parameters and axes parallel to x = o will 
pass respectively through p and q when their vertices are at 
N and M, the feet of the ordinates at q and p. It then becomes 
evident that the figure will be greatly simplified by a trans¬ 
formation which corresponds to the substitution 

x = X ; y\ = Y 

i.e. by raising each ordinate to the power -J. The derived curve 
ADB (fig. 1) which replaces the cycloid has for its equations 

X = 0 — sin 6 y Y = (1 — cos Of = 2 \/2 sin 3 ... (4) 

The semi-cubical parabolas are replaced by the straight lines 
PN, MQ, inclined to AB at angles a, ft, such that 

tan a = Gp>j2s/2kt\ tan ft = <7^/2 \/ 2kt ... (5) 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjoumals.org/ at University of Massachusetts Medical School on April 4, 2015 










dLV f’£9 '•syhnw£06i 


Jan. 1903. . Mr. Aldis’Tables of i(d+cos 0). -| I49 

4. A graphical method of finding the required approximations 
is thus suggested. The curve ADB given by (4) must first be 
constructed on squared paper, and the base line AB may 
conveniently be graduated in 0 . Then it is necessary to place 
the lines PN, QM at known inclinations to the base in such a 
manner that their intersections with the curve project (as in 
fig. 1) on to their intersections with the base. A convenient 
way of doing this is to draw the lines PN, QM on ground <dass 
with a third line in the direction of the base, and to move the 
figure so obtained over the paper until as good a fit as possible is 
obtained. Bough values of e and S can then be read off. But 
actual trial of the method in one or two cases has seemed to 
show that it is not susceptible of any great accuracy. This is in 
a great measure due to the fact that the given interval of time 
is in general small. The analytical method of approximation 
now to be described takes advantage of this fact, and will 
become the more efleetive as the graphical method becomes more 
impracticable. _ . 

5. In fig. 2, which Corresponds to the essential part of fig. j, 



™- PK b , e drawn Parallel to MN, and let the angle QPK = v 
ihis angle is easily determined, for r 


tan y = tan /3 — tan a 

Let T be the point of contact of that tangent to the curve which 
that ^ d *° and l6t G = 7 at this P oint ‘ Ifc is easily found 


tan y = 


- 4 = cos \ T 

v 2 
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150 Mr. Plummer , Note on the Use of lxiii. 3, 

Hence by (5) 

cos i T = (crfi — crfi)/Gkt . (6) 

a result which presents a curious analogy to Euler’s theorem. 
In fact, the latter is hereby deduced from Lambert’s theorem * 
if it is noticed that, in the limiting case of the parabola, 

6. It will be convenient to pub 

2/? = f + S; 2^ = £— 8 . (7) 

Now r must clearly lie between 8 and e, and is indeed a good 
first aproximation to the value of //. Let L be the corre¬ 
sponding value of £ which makes (2) and (3) consistent. Bub 
these equations give 

(.Tj sin \ S = o' 2 sin ± e 
which can be written in the modified form 


tan U = 


*1 + 01 


tan ^7 


Hence 


tan \£ x = 1 —- tan \r 
i+n 


( 8 ) 


( 9 ) 


where n = o^/ov Thus r and £ x are found very easily by (6) and 
(9), and the corresponding approximations to 8 and s are deduced 
immediately. 

7. It is now easy to see how the approximation can be im¬ 
proved by a simple step. Let (fig. 2) i be the foot of the ordi¬ 
nate through T, m the foot of the ordinate midway between PM 
and QN, and h the foot of the ordinate through the point for 
which 0 = \ (e + 8). Let y = <f> (x) be the equation of the curve 
and x and x -j- h the abscissae of P and Q. Then by a well-known 
theorem 

wtiprf* /* — X -J- X ^ 7. 

where j _ 2 + 24 . A 


if the expansion of / is not carried beyond the first power of x. 
Hence 


tm = 


24 


f"( x ) h -z 


* Adams has given a simple proof that Euler’s theorem cm be deduced as 
a limiting case.— Collected Scientific Papers, vol. i. p. 413. 
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And if x = ip(d) and 6 . becomes d+Af) at Q, then 
km — |^(d)+|-i//(d + A6>) — ^(d+±\ 6 ) 

= iV'(6) . \e> 

••• At = {if'(«)+* • ^ tmy } AS* 

= + sin 0 . A 0 2 

when the differentiations are performed. The correction to be 
applied to r in order to give a second approximation to the value 
of ?/ is therefore 

—ht -4- i//(0) = — cot . A0 2 

If, as above, the calculation is made in terms of the lower end of 
the chord, 0 = S ; and starting from the upper end we should 
find 0 = £. It is therefore natural to adopt 6 — t. Also 
A0 = 2^. Hence it is inferred that if 


Ar = -t 5 2 cot It • £ 2 .(to) 

then t+A r is the second approximation to the value of from 
which improved values of e and 3 can be deduced in the manner 
already explained. All the angles occurring in the foregoing 
formulae and the constant k are supposed to be expressed in 
circular measure. 

8. The result of practical trials is to show that the method is 
extremely effective. In the first instance the example which was 
taken by Mr. Aldis from the Theoria Motus was tried. This 
gave 


o / n 

t = 78 38 21*6 ; 
r= 5 1 4 ' 8 ; 


Oil, O l ,, 

4 X = 3 8 21*6; fr = 81 46 43*2 
£2 = 3 8 4*° j € 2 = 81 41 10*8 


How Mr. Aldis finds for the correct value of €, 8i° 41' n"^, 
and hence for log a , 0*4224385, while the value given by Gauss 
for log a is 0*4224389. The above value of e 2 gives log a 
= 0*4224400. And if Ar is recalculated by means of in the 
place of <f x , the correct value of e is obtained. A second trial 
was made with an artificial example given by Dr. Moulton * in 
which the eccentricity was unusually large. In this case the 
value assigned to a was 2*65. The logarithms of three radii 
vectores are given only to six figures. The value found for f 2 
from the shorter period of thirty days gives a = 2*650063, which 
shows a fair agreement with the value actually assigned. But if 
the extreme places of the planet, with an interval of sixty-five 
days, are taken, the case is by no means so good. From the 
value found for e 2 , a — 2*6519. It is true that an improvement 
can be made by recalculating Ar as mentioned above, but the 


* Astr. Journal, No. 510. 
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method would lose the advantage of its simplicity thereby. It 
is better to treat as an approximation which facilitates the use 
of Mr. Aldis’ tables. Cases in which the interval of time is so 
long must occur very rarely, if at all, in the practical applica¬ 
tions of the problem. But in such cases, and whenever the 
highest degree of accuracy is required, the value of the tables, as 
far as this problem is concerned, will consist in giving the final 
improvement, or at least the verification, to approximate solutions 
found by some preliminary method. 

9. It is, however, desirable to' examine the nature of the 
second approximation more closely. The tangent to the curve 
(4) at the point 6 =t is parallel to the chord joining the points 
6 —rj—£ and 0 = 7} + £. Therefore 


i cos J- 


, 3 cos \t 


1 sin 3 j(y + 4 )—sins -§ — 

2 . 1 —sin (?; + £; +sin (r/— '() 

3 cos J47 sin -|£—cos -f rj sin § £ 
X —cos i] sin X 


3 cos \r 

2 COS ^77 


-K(3c I -3«3)-6.8^(3 c t-27fl3) + i 20 .^ 3 (3 c i- 2 43«3)-- • • 


£ ( c i c 3 ) +g-<; 3 (ci + c 3 ) 1 1 o 2 5 (^i 4 * £3) + • • • 


where 


C T — COS ^-77 ; C 3 = COS ±7} 


The expansion in X is clearly an even function, and if terms of 
the order £ 6 be neglected 


cos 


-It - c 3 ) -6.4' 2 ( Cl ~9 c s) + iio. Tb^ - 8IC 3) 


c °s In ( Cl —c 3 ) + K 2 (c I + c 3 ) —rbf 4 ( <; 1 + ° 3 ) 

. cos£7—cos -|t_ 2 S jC( c i ““a) - T sW' 4 ( 1 7^—655) 


cos jr 


( Cl —C 3 )—9c,) 


_ s_;'a_ i_ Ci + 2 S 5 c 3 

24^ ‘ Vs 

6 1 °3 

= 21 f 2 + J fVo ( 1 9 1 - 64 cot 2 ■§!/)£+ 

10. In the coefficient of r may be written for 77. Then 

sec Jr A cos Jr = srW^ 1 “ 6 4 cot 2 Jt);4 

where A cos Jr is the correction to be added to cos Jr to obtain 
cos J 77, so that 

A cos 1-7 = —J sin Jr. At — -J cos -Jr(Ar) 2 
Hence, to the degree of approximation adopted, 

A- = —T^cot tr. i 2 -2 ? Vo c0t i T ( [ 9 I + 61 COt 2 -|r),4 
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which confirms the formula (10) previously obtained and gives 
a further development. Now since 

4 tan 2 K=£ 2 + K 4 

it is possible to write 

Ar=-f cot |r tan 2 ■§■£- cot M 61 cot2 i' — 9>f 4 ■•• (11) 

Therefore if powers of f above the third be neglected it is more 
accurate to write (10) in the form 

Ar=—■§• cot \t tan 2 £ 

which again, by means of (9), may be expressed 

AT =-i(f^Y Un ? T . < I2 > 

Thus the second approximation to y may be derived directly 
from (6) and (12) before an approximation to f has been calcu¬ 
lated. The work to this stage will therefore be very rapid, and 
it will be easy to judge from the value then found for £ whether 
it is worth while to refine on the approximation by means of 

(”)• 

11. The representation of the consequences of Lambert’s 
theorem by means of a curve related to the cycloid was suggested 
by the method I have previously described * for finding graphi¬ 
cally an approximation to the solution of Kepler’s equation by 
means of the trochoid. The application of Mr. Aldis’ tables to 
the latter problem will be rendered more advantageous by using 
them as an auxiliary to some such approximate method as the 
one to which I refer. My description of that method was not, as 
it soon appeared, original. It had, for instance, been published 
by Dr. Bambaut + a few years earlier; and quite lately it has 
again been brought forward by B. Gonggrijp. X But the 
principle of the method is as old as Wallis and Newton. 


Note added 1903 Jan. 24. 

The purpose of the foregoing note was, as explained in the 
beginning, merely to supply an approximate method which might 
serve to facilitate the use of Mr. Aldis’ tables, and the 
examples mentioned in § 8 suffice to show that this object has 
been fully attained. But the method as further developed in 
§§ 9 and 10 is something more than a mere auxiliary to the 
tables. It provides an adequate and convenient means of solving 
Lambert’s equation in cases of moderate difficulty without the 
aid of any special tables. 

The essence of the method is contained in the equations (6), 

* Monthly Notices , vol. lvi. (1896), p. 317. 
t Ibid. vol. 1. (1890), p. 301. 

T ! X A. N. (1901), No. 3720. 

M 2 
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(8) and (11), by means of which the values of 77 and 4 are found 
by successive approximations. It is well to remark that 

— E) ; cos 77 = e cos \ (E' +E) 

where E, E' are the eccentric anomalies of the planet; and that 

ap — rr ! sin 2 A cosec 2 £ 

wffiere 2 A is the difference of the true anomalies. This explains 
the nature of the expansion (11), which might, if necessary, be 
developed further without serious difficulty. Its advantage con¬ 
sists in its rapid convergence for moderate values of £ and in the 
fact that to the order & it consists of only two terms, of which 
the second is very small. For actual calculation (11) may be 
written in the form 

—iAr = Z + FZ 2 . (iE) 

where Z = cot \r tan 2 -J-f 

F = zini 61 —9 t an “ i T ) cot ±r 

When log Z is found, Z can be taken from a table of logarithmic 
sines ; and the following table, which will apply to eccentricities 
of at least 0*5, gives log F in such a way that the term FZ 2 will 
be found in seconds of arc :— 





Table of log F. 




It- 

log F. 

Diff. 


log F. 

Diff. 

It. 

log F. 

Diff. 

O 

30 

4-9186 

— 192 

O 

40 

473°S 

— 190 

O 

50 

4'5238 

— 240 

31 

•8994 

190 

4i 

7H5 

193 

51 

•4998 

250 

32 

•8804 

189 

2 

•6922 

195 

52 

•4748 

262 

33 

•8615 

187 

43 

•6727 

198 

53 

•4486 

276 

34 

•8428 

187 

44 

•6J29 

202 

54 

•4210 

29s 

35 

•8241 

186 

45 

•6327 

206 

55 

•3919 

310 

3*3 

•8055 

187 

46 

■6121 

212 

56 

•3609 

332 

37 

'7868 

186 

47 

•5909 

217 

57 

•3277 

358 

3* 

•7682 

1S8 

48 

•5«92 

223 

58 

•2919 

390 

39 

47494 

- 189 

49 

4-5469 

-231 

59 

•2529 

-430 







60 

4-2099 



The more complete formulae have been applied to the example 
already mentioned in § 8, in which the interval of time is 65 days 
and the eccentricity 0-3. The result is quite satisfactory, and 
shows that the formulae are as accurate as the data require. 
As the second interval between the observations does not enter 
into this part of the computation, and is therefore unlimited, 
this example is a far more severe test than when two intervals 
of 30 and 35 days are used. The additional term FZ 2 amounts 
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to 4"‘7 7. It may be necessary, as it is in this case, to repeat 
the calculation by means of (8) and (11') two or three times, 
but the process is in practice very simple and rapid. It is 
convenient to use logarithmic tables giving every second of arc 
(e.g. Bagay’s), to keep the openings for i/7, if, and iAr marked, 
and instead of repeating the whole of the figures (which in any 
case are very few) it is sufficient merely to write the corrections 
in parallel columns. In fact, when rr x and /r 2 are known it takes 
very little time to find 77 and £, and hence a, p and e. The 
method is summarised in the following scheme :— 


Summary of Calculation . 

(1) Calculate tr„ o* 2 by means of 

tr 1 * = r+r , + c ; <r 2 2 = r-fV — c 

(2) Calculate w = o , 2 /<r I , and hence log N = log (1— n) 

-log(i +n) 

(3) Find \t by means of 

cos \r = (r t-P —<t 2 3 ) / 6kt 

(4) Find iAr from 

iAr=-^tanir 

(5) Hence i/7 2 = i(r +Ar), and f 2 can be found from 

tan if 2 = N tan i/7 2 

(6) With f 2 as an approximate value of f find log Z from 

^ = 4 cot ir tan 2 if 
and recalculate iAr by means of 

—iAr = Z + FZ 2 

(7) If necessary repeat the operations (5) and (6) until the 
equations 

* —i(/7—V) = —iAr = Z+FZ 2 

tan if = 1 ST tan i?7 

are both satisfied. 

(8) When 77 and f are thus known with sufficient accuracy 
a is found from 

a = <r J 2 / 4sin 2 i(i/ + £) 

and p from 

ap — rr } sin 2 A cosec 2 £ 

Finally the eccentricity is found from a and p. 

University Observatory , Oxford : 

1903 Jan. 8. 
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On a New and Accurate Method of Determining fame, latitude, 

and Azimuth with a Theodolite. By W. BT Cooke, M.A., 

Government Astronomer, Western Australia. 

The principle of this method is not new. It was, I believe, 
first advocated by Mr. S. C. Chandler, who designed an instru¬ 
ment which he called an “ Almucantar, ,; developed the mathe¬ 
matical theory, and made a series of very fine observations at 
the Harvard College Observatory. An instrument of that kind, 
however, does not form part of a surveyor’s ordinary outfit, 
whilst a theodolite does ; and my present purpose is to show how 
the method may be applied to this universal instrument, and how 
by its means results of far greater accuracy than those usually 
obtained may be deduced. Briefly, the spirit level can be used 
instead of a mercury flotation in order to ensure constancy in 
zenith distance, and this brings the method within the power of 
every surveyor. 

I shall first show it in connection with a 5-inch theodolite, 
and state the results of a few observations already made. After¬ 
wards I propose to show that with a 12-inch instrument results 
can be obtained which, I believe, exceed in accuracy anything in 
field work hitherto published. 

Fire-inch Theodolite . 

The observation for time and latitude consists in observing 
the time of transit of Nautical Almanac stars over the hori¬ 
zontal thread when the instrument is set so as to sweep a small 
circle in the sky parallel with the horizon, and at an altitude 
about equal to the observer’s latitude. For azimuth it will be 
necessary, in addition, to take the time of transit across the 
vertical thread and to read the azimuth circle. An ordinary 
watch which possesses a seconds hand will do for timekeeper. 
The advantages may as well be enumerated here. 

For time and latitude errors of vernier reading do not exist, 
as we require the circles for approximate setting purposes only. 
Extreme accuracy of construction and adjustment is unnecessary. 
We require only ordinary care in levelling. Error of collimation 
is absolutely immaterial. We do not even require that the 
bubbles shall be properly adjusted. All we require in the way 
of adjustment is that the vertical axis shall be reasonably ver¬ 
tical, so that the cross bubble on the altitude circle shall remain 
fairly steady as the instrument is swung round in azimuth, and 
great accuracy is not required even for this. As a matter of 
fact the instrument I used had been standing in a corner for 
months, and I simply took it out and levelled it up. 

Yet, notwithstanding this apparent want of respect, results 
of considerable accuracy may be easily obtained. Those obtained 
in two evenings’ trial will be stated almost immediately. 
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